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On guided electromagnetic waves in photonic 
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Abstract. The paper addresses the issue of existence and confinement 
of electromagnetic modes guided by Unear defects in photonic crystals. 
Sufficient condition are provided for existence of such waves near a 
given spectral location. Confinement to the guide is achieved due to a 
photonic band gap in the bulk dielectric medium. 

1. Introduction 

A photonic crystal, also called photonic band-gap (PBG) material, is 
a periodic medium which plays the role of an optical analog of a semi- 
conductor. Such a medium has a gap in the frequency spectrum of electro- 
magnetic (EM) waves. The idea of a photonic crystal was first suggested 
in 1987 [16,30], and has since been intensively studied experimentally and 
theoretically (see, e.g., the recent books [15, 17,25,26,31], the mathemati- 
cal survey [20], the on-line bibliography [23], and references therein). This 
interest has been triggered by the numerous promising applications of PBG 
materials, one of which is using photonic crystal for manufacturing highly 
efficient optical waveguides. The idea is to introduce a linear "defect" into 
a PBG material, and to guide through it EM waves of a frequency prohib- 
ited in the bulk. Numerical and experimental studies have shown that such 
superior guides can be efficiently created, e.g. [15, 17,22,25,26]. 
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In order to create such a guide, one needs to establish several facts. 
The first, and foremost, is existence of guided waves of frequencies in the 
band gap. The second, and an easier one, is confinement of these modes. 
This paper addresses both issues, by finding some sufficient conditions of 
existence of guided modes, and showing their confinement to the guide, in 
the sense of being evanescent in the balk. Similar results were previously 
obtained by the authors in [21] for scalar models (i.e., for acoustic analogs 
of PBG waveguides). In this paper, we will address the above questions for 
the full Maxwell case. 

There is another important question to be resolved. Namely, one needs 
to show that the impurity spectrum that arises in the spectral gaps due to the 
presence of a linear defect does not correspond to bound states. This diffi- 
cult issue is not addressed here (see some relevant remarks and references 
in Section 5). 

In Section 2 we introduce the main model to be investigated. The next 

Section 3 contains formulation of the main results, with the proofs provided 
in Section 4. The paper ends with the sections devoted to final Remarks and 
Acknowledgments, as well as the Bibliography. 

2. Preliminaries 

We start by describing the mathematical model studied in this paper. 
Let eo{x) be a bounded positive measurable functions in separated from 
zero: 



(2.1) < Co < £0(2;) < ci < 00. 

It is usually assumed in photonic crystal theory that £0 is periodic with 
respect to a lattice F C M^, but this is not required for our results. 

The function eq represents the dielectric properties of the bulk material. 
In other words, one can think of the space filled with a dielectric material 
with the dielectric function eq. 

The unperturbed Maxwell operator Mq is the self-adjoint realization of 

(2.2) Mo:=V^^-V^ 

Eo{x) 

in -L2(M^; C^) defined by means of its quadratic form 

(2.3) j EQ^\Vu\^dx 

with the domain H^{M?] C^). We use here the shorthand notation 

V^u = V X 14 = curl u. 
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A cylindrical domain Si (see Fig. 1) will represent a linear "defect 
strip": 

Si := {x = (xi, x') eR^\xieR,x' e in}. 
Here the cross-section m of the strip is a domain in (e.g., the unit ball 
centered at the origin), scaled with factor /. 




Figure 1 . The waveguide Si . 



(2.4) e{x) 



We now introduce the perturbed medium with homogeneous dielectric 
properties inside the defect strip Si. 

£ > for X G iS; 
(x) for X ^ Si 

The perturbed Maxwell operator 

(2.5) M:=V''^V^ 

e{x) 

corresponds to the medium with the linear defect. It is defined, analogously 
to Mo, as a self-adjoint operator in L2(]R'^; C'^). 

Remark 2.1. The reader has probably noticed that we disregard the 
standard restriction of the Maxwell operator to divergence- free fields [15] . 
The difference is essentially in acquiring a huge eigenspace corresponding 
to the zero frequency, with all other parts of the spectral decomposition 
staying intact. Since the problem of guided waves concerns the situation 
inside the spectral gaps of the operator Mq, this difference is irrelevant in 
this case. On the other hand, abandoning the zero divergence condition will 
simplify the techniques considerably. 

Our goal is the same as in the paper [21] devoted to the scalar (acous- 
tic) case, i.e. to show that for any gap (a,/?) in the spectrum cr(Afo) of 
the unperturbed medium, under appropriate conditions on the parameters 
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/ and e of the line defect, additional spectrum arises inside the gap, with 
the corresponding (generalized) eigenmodes being confined to the defect 
(evanescent in the bulk). 



In order to formulate our first main result, we need to introduce the 
following quantity: 

Definition 3.1. We denote by > the lowest eigenvalue of the 
Laplace operator A acting on divergence free R^-valued vector fields on Q, 
with Dirichlet boundary conditions on dO,. 

I.e., z/ > is the smallest number for each a non-trivial solution of the 
following problem exists: 



where ^ is the external normal derivative on 30,. 

Our main results are given in the following theorems. 

Theorem 3.2. Let G = {a, (3) be a non-empty finite interval, such that 
a > Q (we will be especially interested in the case when G is a gap in the 
spectrum of the "background medium" operator Mq). Let the following 
inequality be satisfied: 



Then the interval G contains at least one point of the spectrum a{M) of the 
perturbed operator 

This theorem guarantees that when (3.2) is satisfied, eigenmodes of the 
perturbed medium do arise in the spectral gaps of the background medium. 
Furthermore, if 5 > is such P6£ > u, the corresponding spectrum forms 
a 5-net in the gap. 

Before one can fully associate these modes with the guided waves, one 
needs to establish their confinement to the waveguide (i.e., their evanescent 
nature in the bulk of the material). In order to describe the corresponding 
result, we need to remind the reader some notions and results about gener- 
alized eigenfunction expansions. Here by generalized eigenfunctions one 
understands solutions of the eigenvalue problem, which do not decay suffi- 
ciently fast (or do not decay at all) to belong to the ambient Hilbert space L^. 
One can find detailed discussions of generalized eigenfunction expansions, 
for instance, in [3, 4]. 



3. Formulation of the results 



(3.1) 




(3.2) 
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Namely, as it is shown in [18], for the operator M that we consider, for 
almost any (with respect to the spectral measure) A in the spectrum a{M), 
there is a generalized eigenfunction ux{x) e Hl^^{W?, C^) with the growth 
estimates 
(3.3) 

(l + \x\)-^u{x)eL2{R^-X^) and (1 + V^m(x) e L2(M^; C^*) 

for some iV > 0. This system of generalized functions is complete in the 
whole space. For elliptic operators with smooth coefficients, this is a well 
known fact [3]. 

Definition 3.3. We refer to the following growth condition as poly- 
nomial boundedness of order N: for any compact set X C and x e R^, 

(3.4) \\u\\l2{{K+x);C^) + II V'''u||L2((i^+^);C3) < Ck{1 + 

In the next result we will use the following notation: for x = [xi, x') e 
R^, we denote by Xx{y) the characteristic function of the cube 

{y I \yj-Xi\ < iforj = 1,2,3} 

centered at x. I.e., Xx{y) is equal to 1 when y is in this cube, and other- 
wise. 

Theorem 3.4. Let G be a finite spectral gap of Mq and u\ be a poly- 
nomially bounded generalized eigenfunction of M corresponding to X & 
G n a{M). Then there exist positive constants G andG{X) such that 

(3.5) llxxt^All < C(l + |xi|)^e-^W'^-*(-'^'), 
where N is the order of polynomial boundedness ofu\. 

When the bulk medium is periodic in the -direction, the polynomial 
growth in (3.5) disappears: 

Theorem 3.5. If eo{x) is periodic in the xi-direction, then one can 
find a complete family of generalized eigenfunctions that satisfies 

(3.6) \\u\\\l2{{k+x);C-') + \\^^ux\\l2{{k+x);C^) < Gk{1 + 

for any compact set K C R^, x G R'^. In this case, for A e G n cr(M), one 
has the estimate 

(3.7) ||x.«A||<^^e-^W^'*^(^'^'). 

4. Proofs of the results 

In what follows, the norm and inner product in L2(R"; C^) will be de- 
noted by II • II and (•, •) respectively. 
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4.1. Proof of Theorem 3.2. We will show that if > and 5 > is 
such that PSe > u, there is spectrum of the operator M in the (5-vicinity of 
fx. Then, taking fj, — {a + P)/2 and 5 — {P — a)/2, one gets the statement 
of the theorem. 

In order to show this, due to self-adjointness of M, it is sufficient to find 
a vector function w e -L2(M^, C^) of unit norm, such that 

(4.1) \\Mw - i2wf< 5^. 

Let 5f be a smooth, unit L^-norm, divergence free real vector field on 
with compact support in n and unit -norm, i.e, 51(1/, 2;) = z), -2)) 
where (f),Ce C^{n) and V • ^ = 0. We define 

gr.= {<j)iXi) = l-\H^'/l)X{x'/l)). 

Then gi is also a unit L^-norm divergence free field. Let ip{xi) e C^(R) 
have unit L2(M)-norm and il^nixi) = n~^^'^ip{xi/n) for n > 0. Clearly 
tl^n{xi) also has unit L2-norm. 

Denoting k = ^JJie, we introduce the following candidate for an ap- 
proximate eigenfunction: 

{ ' 

(4.2) wi,n{x) = M^iy""' Mx') 

\ CM 

The function wi^n clearly has unit norm in L2(M x l^l). 

Instead of estimating the left hand side of (4. 1), we will estimate 
II e {Mw — jjLw) I p. Taking into account that the function w is supported 
inside the defect, the needed inequality (4.1) can be also rewritten as 

(4.3) ||V^V^«;-A;'HI'< 
Using the identity 

V^V^w = -Aw + V(V -w) 
and that gi is divergence free, we obtain 

||V^V^w-A;2w| 

where the norms are in L2(M^; C^). 

Since the functions 0, C and i}} are real valued, their assumed normal- 
ization shows that the above expression is equal to 

n-' irilL(M) +^k'n-' II^IIL(M) l|A<7|lL(a,M.) 
+2(n0 '0)l2(m)(A^, g)L^(Q,M?y 






■« + 2^/cO0;-7/'„A0z 

-(C + 2z^O0-^nA0 
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Since n can be chosen arbitrarily large, the terms with the factors that are 
negative powers of n can be made arbitrarily small (uniformly with respect 
to k on any finite interval). Hence, one needs to control only the remaining 
terms by an appropriate choice of a divergence free vector field g. In other 
words, one is interested in making l^'^ WAgWj^^^^ smaller than i.e. 

(4.4) ||A^|U,(0;R2)< l^Ss, 

while keeping {Ag, g) under control. 
Since 

V = inf ||Ac/||i2(n;R2), 
where the infimum is taken over real, unit L^-norm, divergence free vector 
fields g e C^(Q; M^), this condition boils down to 

(4.5) f 5£ > V, 

which proves the statement of the theorem. □ 

4.2. Proof of Theorem 3.4. We assume here that A belongs to a finite 
gap G of the spectrum of Mq and it := it a is the corresponding generalized 
eigenfunction of M satisfying (3.4). Let us also introduce the resolvent 
R{\) = (Mo — A)~^. We will also use the function Xx introduced before 
Theorem 3.4. 

We will need the following auxiliary statement concerning the exponen- 
tial decay of the resolvent, which is a result of [8]: 

Lemma 4.1. [8] There exists a positive number m\ that depends only 
on the distance of the point \from the gap edges, such that for a positive 
constant C, the following estimates hold for the local L2(M^; C^)-norm of 
the resolvent R{X): 

^ ^ ^ ||x.VXi?(A)x.||<Ce— 

for any u, v E Ml Here the norms in the left hand side are the operator 
norms in -^2(1^^; C"^)- 

We now consider the sesqui-linear form 

Q[ip,w] := (VV, -V^i/;) - X{ip,w) 

with the domain H'^{R^; C^). 

Let (p :— R{X)xxU. Note that (p belongs to the domain of the operator 
Mo. 

Letp = max (2dist(a;, Si),l) and ^x{y) be a nonnegative smooth cutoff 
function that depends on yi only, is supported in (xi — (p+ 1), xi + (p+ 1)) 
and such that it is equal to 1 on [xi — p,xi + p]. We assume further that 
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ix{y) < 1 and \V^x{y)\ < C for some constant C and all x,y e M^. Note 
that ^^u e if^(R^; C^). Using w = ^^u, one gets 

Thus, the goal is to estimate Qi^?, Cxu] from above. On the other hand, using 

the equality Mu — Xu and integration by parts, one obtains 

(4.7) 

= (V>,eVx(e.H)) + {V^^,e-'V^Cxu)) - {V^Cx^),e-'V>'u), 
where we used the notation 

.^j, L 

eo{x) e{x)' 

Notice that e is supported inside the strip Si. 

Using the identity V^{^u) = ^V^u + x u, the first two terms in the 
last line of (4.7) can be combined to obtain 

(Vx<^, e(e.VXM + VCx X u)) + (V>, e-'iCxV^'u + V^. x u)) 
= (VXv9,e.eVXM) + (V>, (e + 5-i)Ve. x u) + {V'<ip,e-'^,V'<u) 

The term —{V^{$,x'^),£~^V^u) can be expanded to 
Combining the last two expressions, we get 

(4.8) (VX(^,e.eVXH) + (VX(^,£-ive. xu)- (V^. x ip,e-^V''u) 

Our last task in proving the theorem is to estimate from above the terms 
in (4.8). Let V = [xi — p— 1. xi +p+ 1] x ICl. This is a compact domain that 
can be covered by the union of p fixed size domains Vj = [aj , aj+2] x IQ and 
which contains the supports of {^e). Also note that dist(a;, Vj) > dist(a;, Si). 
Using Lemma 4.1 and (3.4), we get for any < 77 < nix 

^, <c ||E,xv.vxit:(A)x.«||||E,x^.vx«|| 

(4.y) <CpWxi\+p+ iyN^-rn,dmx,S,) 

< C{\xi\ + l)2^e"("*^"'')^^^^(^''^'^ 

We used here that p — max(2dist(x, Si), 1) and denoted by C different 
constants. 

Let us move now to estimating the last term in (4.8). Denote by a > 
a number such that shifts of IQ by vectors aj with j e cover the whole 
space R^. We denote 

Wj := ([xi -p-l,Xi-p\U [xi + p,xi+p + l]) X (IQ + aj) . 
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Then Wj = Wq + (0,aj). Notice that W = UjWj covers suppV^ and 
dist(x,Wj)>C,(p+\j\)-C2. 

We are now ready to estimate the last term of (4.8) from above. We pro- 
ceed as before, using the lemma, the polynomial growth of u, and uniform 
boundedness of V^x- 

j 

<CJ2i\xi\+p+\j\ + l)2iVe-m,dist(a.,W-) 

(4.10) j 

< Ci\x,\ + i)2iVg-c;m,dWA), 

The middle term in (4.8) is estimated analogously. Combining these 
estimates, we get 

\\XM\'- Q[V,^.U] < C{1 + |^,|)2iVg-C.n.dist(.,50_ 

This finishes the proof of the theorem. □ 

4.3. Proof of Theorem 3.5. In this periodic situation, operator M has 

a complete family of generalized eigenfunctions that do not grow in the 
xi-direction of periodicity. Indeed, according to Bloch-Floquet theory [19, 
24], a generalized eigenfunction u = ux of M corresponding to A can be 
chosen as u{x)e^''^^, where u{x) is periodic in xi-direction with period a 
and k — 27r/a. Thus, u satisfies (3.6). Then, repeating the previous proof, 
one comes up with the estimate (3.7). □ 

5. Remarks 

(1) Theorem 3.2 provides sufficient conditions for the existence of a 
5-net of the defect spectrum inside a spectral gap. One wonders 
how much of the gap the guided mode spectrum can occupy. Can 
it fill the whole gap ? If so, under what conditions ? There seem to 
be no rigorous results available concerning these questions. 

(2) Results of [2] on improved Combes-Thomas resolvent estimates 
show that the exponential decay constant mx, which clearly de- 
pends on the distance of the point A from the spectrum, behaves as 
i/(A — q;)(/5 — A) inside a gap G = {a, j3). 

(3) As it has already being mentioned, in order to have the full right to 
call the discovered modes "guided", one needs to show that they 
do not correspond to point spectrum (i.e., to bound states). Here 
the most treatable case should be of a periodic medium with a lin- 
ear defect aligned along one of the lattice vectors. In this situation 
one can apply the Floquet-Bloch theory with respect to the axial 



10 



PETER KUCHMENT AND BENG-SEONG ONG 



variable of the waveguide and hope to use standard techniques ap- 
plied in the case of Schrodinger operators with periodic potential 
(e.g., [5,13,19,24]). This happens to be not an easy task. Even 
in the case of "hard wall" periodic waveguides, when waves are 
contained in a periodic waveguide by Dirichlet, Neumann, or more 
general boundary conditions, this problem is non-trivial. Although 
it has been considered for rather long time [6, 19], the first real ad- 
vances are very recent [11,12,14,27-29]. The case of photonic 
crystall waveguides is more complex, due to absence of complete 
confinement of the waves, which exponentially decay into the bulk, 
but do not vanish completely. Apparently, the most recent work de- 
voted to this issue assumes the bulk material to be homogeneous, 
with an embedded periodic guide [10]. 
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